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Abs  tract. 

The  higher  modal  dynamic  plastic  response  of  impulsively 
loaded,  fully  clamped  beams  is  examined  herein,  using  various 
rigid  perfectly  plastic  theoretical  procedures  and  a numerical 
elastic-plastic  computer  code. 


-ii- 


Hotation 

p pressure 

t time 

t*  defined  by  equation  (47c) . 

u axial  displacement 

w lateral  or  transverse  displacement  defined 

in  Figure  1. 

x,z  axial  and  transverse  coordinates  defined  in 

Figures  1 and  3,  respectively 
B,H  beam  breadth  and  beam  thickness 

Er  energy  ratio 

2L  beam  length 

M,N  bending  moment  per  unit  length  and  axial  force 
per  unit  length  defined  in  Figure  3. 

Mq  o„H2/4 

N OH 

o • 

Q transverse  shear  force  per  unit  length  defined 

in  Figure  1. 

q o H//T* 

O 0 

Vo  peak  value  of  initial  impulsive  velocity 

W transverse  displacement  at  an  interior  hinge 

Hp  final  or  permanent  transverse  displacement 

W#  If/H 

W*f  Wf/H 

X defined  by  equation  (21b) 

I*  mass  per  unit  length 

p density  of  materiel 


-iii- 


o.l%  uniaxial  static  proof  stress 
3 ( )/  3 t 

3 < )/  3 * 


1.  Introduction. 


The  higher  modal  dynamic  plastic  response  of  fully 
clamped  beams  subjected  to  large  impulsive  loads  was 
examined  recently  as  described  in  Reference  {1 ] . Exact 
theoretical  solutions  were  presented  for  the  first  three 
modal  responses  of  rigid  perfectly  plastic  beams.  It  was 
evident  from  some  experimental  tests  conducted  on  aluminum 
6061T6511  beams  and  reported  in  Reference  [1]  that  geometry 
changes,  or  finite-deflections,  played  an  important  role 
in  the  dynamic  response,  which  agrees  with  previous 
observations  on  the  behavior  of  uniformly  loaded  beams 
(e.g.  [2]).  Thus,  the  theoretical  procedure  developed  in 
Reference  1 3 ] was  used  to  examine  the  influence  of  finite- 
displacements  on  the  fundamental  modal  response. 

The  infinitesimal  and  finite-deflection  analyses 
developed  in  Reference  {1]  are  further  generalised  in  order 
to  predict  the  behavior  of  impulsively  loaded, fully  clamped, 
rigid  perfectly  plastic  beams  which  undergo  any  symmetric 
or  antisymmetric  modal  response.  The  numerical  elastic- 
plastic  behavior  of  fully  clamped  beams  subjected  to  impulsive 
"modal"  velocity  fields  is  also  examined,  using  the  spatial 
finite-element  JET  3C  computer  program  of  Wu  and  Witmer  (4,  5] 
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2.  Fundamental  Equations. 

The  governing  equation  for  a rigid  perfectly 
plastic  beam  which  undergoes  a modal  response  of  the  form 

w(x,  t)  - ♦(x)W(t)  (1) 

is  M"  = -u$W  (2) 

provided |m )<_  Mq  and  where  the  various  symbols  are  defined 
in  the  accompanying  Notation  and  Figure  1.  The  boundary 
conditions  at  the  fully  clamped  supports  (x  = 0)  of  a beam 
are 

$ = 0 , M ■ +M0 
and  0=+  1,  M = +M0  , Q = 0 

at  x a L of  a beam  undergoing  a symmetric  mode,  while 

* « 0 , M - 0 (5a,  b) 

at  the  centre  (x  » L)  of  a beam  responding  in  an  antisymmetric 
modal  form. 

3.1  Symmetric  Modal  Response  (Infinitesimal  Displacements). 

Any  symmetrical  modal  response  t(x)  * 4>  (2L  - x) , which 
satisfies  equations  (3a)  and  (4a),  consists  of  two  types  of  shape 


(3a,  b) 
(4a-c) 


functions 


3 


<J>,  = +x/n,  , o < x < n. 


in  region  1 adjacent  to  the  fully  clamped  supports  shown  in 
Figure  2(a)  and 

_ (n^  - x)W.  . (x  - n.  ,) 

*i  = + <it  - t (n,  - ^ ' "i-i  i * i "i  (7) 

for  any  region  i,  where  2 <_  i £ n and  2n  - 1 is  the  mode 
number.  The  upper  signs  - and  + in  equation  (7)  refer  to 
any  region  similar  to  region  3 in  Figure  2(a) , while  the  lower 
signs  + and  - are  associated  with  any  region  similar  to 


region  2. 


Now,  substituting  equation  (6)  into  equation  (2)  gives 


M"  = +yw1x/n1 


provided  0 <_  x Hj  . It  may  be  shown  by  integrating 
equation  (8)  twice  with  respect  to  x that 


W,  « -eMo/fur^2) 


when  satisfying  equation  (3b)  for  the  plastic  hinge  at 
x • 0 and  using  the  requirements  M = +Mq  and  Q = 0 for  the 
plastic  hinge  at  x ■ n a to  evaluate  the  two  constants  of 
integration. 
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Similarly,  substituting  equation  (7)  into  equation 
(2),  then  integrating  once  with  respect  to  x and  satisfying 
0 ■ 0 at  the  plastic  hinges  located  at  x = and  x = 

leads  to  the  requirement 

• • 

wi-i  “ V (10) 

A further  integration  of  equation  (2)  with  respect  to  x 
gives 

(T1i  ' ni-l)2  * -12Mo / (nWi)  (11) 

when  using  equation  (10)  and  the  conditions  M = +Mq  at 

x = ^i-i  an<*  M “ +Mq  at  x = ni,  where  the  upper  and  lower 

signs  of  Mq  are  associated  with  the  upper  and  lower  signs 

in  equation  (7) , respectively.  Equation  (10)  is  valid  for 

• • 

2 i i 1 n so  that  Wj  = W2  when  i=2  and  therefore  the 
magnitude  of  the  velocity  is  identical  at  every  plastic 
hinge  (1/  2,  3,  . . . , n)  and  is  henceforth  denoted  as  W. 

It  is  now  evident  that  equation  (9)  also  applies  to  every 
hinge.  Thus,  no  subscripts  are  required  for  W. 

Equations  (9)  and  (11)  predict  that 

+ /In,  , i - 2,  3,  . . . , n (12) 


+ /Irij  - 


which  gives 


L 


(13) 
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s ince  n = L when  i = n , and 
n 

ri  = n /?n  (14) 

1 2 i 

when  i = 2 . 

However,  it  is  evident  from  equations  (12)  and  (13) 

that  n ■ + j/3n,  = L,  or 
n-3  i 

ni  + (n  - i) /2r\  ^ = L (15) 

when  n - j = i.  Thus,  equation  (14)  and  equation  (15)  with 
i = 2 predict 

n,  = L/(l  + (n  - l)/2}  (16) 

so  that  equations  (15)  and  (16)  give 


ni  _ 1 -i-  (i  - l)/2 
“T  “ 1 + (n  - li/I  ' 


i = 1#  2^  • • 


n (17) 


It  now  remains  to  solve  equation  (9)  for  a fully 

clamped  beam  of  length  2L  subjected  to  an  impulsive  velocity 

■ 

field  which  has  a shape  of  the  same  form  as  Figure  2(a)  with 
equal  peak  initial  velocities  I vj  at  every  plastic  hinge 
1,  2,  . . . , n.  Thus,  integrating  equation  (9)  twice  with 
respect  to  time  gives 


W = VQt  - SMotVfvnj2) 


(18) 


V 


6 


when  using  W = 0 and  W = VQ  at  t = 0 to  evaluate  the 
constants  of  integration.  The  beam  reaches  a permanently 
deformed  profile  at  t = T when  W = 0.  Thus, 

T = VV^V^Mo)  (19) 

and  the  associated  maximum  permanent  displacement  at  all 
plastic  hinges  1,  2,  . . . , n,  is 

W*f  = XnM/12  (20) 

where  W*f  = Wf/H,  X = viVq2L2/ (MqH)  , and  n x = r^/L  (21a-c) 

with  rij  given  by  equation  (16)  for  the  2n-l  mode.  It  is  evident 
from  equation  (1)  that  the  final  profile  of  the  beam  is 

wf  = $Wf  (22) 

Equations  (2)  and  (9)  indicate  that  M"  has  the  same 
sign  as  <f> . Thus,  M"  > 0 in  all  regions  below  the  x axis  in 
Figure  2(a)  since  0 > 0,  while  M"  = 0 on  the  x axis.  However, 
the  boundary  conditions  M = -Mj,  and  M*  = Q = 0 were  used  in 
the  foregoing  theoretical  analysis  for  all  plastic  hinges 
(1,3,...,  n-1)  which  lie  below  the  x axis,  while  M = Mq 
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was  used  for  the  support  hinge  at  x = O.t  It  is  evident, 
therefore,  that  the  bending  moment  distribution  has  minima 
at  the  plastic  hinges  located  at  x = r).  (i  = 1,  3,  . . . , 
n - 1)  and  M > -MQ  for  all  other  parts  of  a beam  with  <f>  >_  0. 

vi 

Moreover,  it  is  now  clear  that  -M0  £ M <_  MQ  within  region  1 
for  which  0 <_  x <_  n j • 

Similarly,  M"  < 0 for  all  regions  which  lie  above 
the  x axis  in  Figure  2(a) , while  M = Mq  and  M'  = Q = 0 at  all 
the  associated  plastic  hinges  (2,  4,  . . . , n) . Thus,  the 
bending  moment  field  has  maxima  at  the  plastic  hinges 
located  at  x = r|.  (i  = 2,  4,  . . . , n)  and  M < Mq  for  all 
other  parts  of  a beam  with  <J>  <_  0. 

It  is  now  apparent  that  the  bending  moment  distribution 
associated  with  the  foregoing  theoretical  solution  has 
-M0  £ M <_  Mq  and  is  therefore  statically  admissible  since  it 
satisfies  the  required  conditions  at  all  the  plastic  hinges 
and  nowhere  violates  the  yield  criterion  in  those  parts  of  a 
beam  lying  between  the  plastic  hinges.  Furthermore,  the 
solution  is  also  kinematically  admissible,  so  that  it  is 
exact  according  to  the  methods  of  plasticity  developed  for 
structures  which  undergo  infinitesimal  deflections. 

Equations  (19)  and  (20)  with  n = 1 and  n = 2 agree  with  the 

+ 

In  fact,  the  theoretical  analysis  also  caters  for  the 
possibility  that  plastic  hinges  (2,  4,6,  . . . , n)  lie 
below  the  x axis,  in  which  case  M = - at  the  support. 
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H 


first  and  third  mode  solutions  presented  in  Reference  [1]. 

It  is  convenient  at  this  juncture  to  list,  for  future 
reference,  the  total  angular  rotation  rates  (8)  across  the 
individual  plastic  hinges  in  a beam  which  responds  in  a 
symmetric  deformation  mode.  It  is  evident  from  Figure  2(a) 
and  equations  (10)  and  (12)  that 


|e0l  " "/»», 


(23a) 


corresponds  to  the  support  hinge,  while 


I © ! I = (ni  + n2 ) w/(n j (n2  - n,)} 


(23b) 


and  | e ± | = 2(ni+1  - (ni+1  - nA)  (ni  - n^) } (23c) 


for  hinge  number  1 and  a typical  interior  hinge  i,  respectively. 


|en|  = I e i I /2 


(23d) 


for  the  central  hinge. 


3.2  Antisymmetrical  Modal  Response  (Infinitesimal  Displacements), 


Three  types  of  shape  functions  <t>  (x)  = -$(2L  - x)  are 
required  to  characterise  any  antisymmetrical  modal  response 
of  a beam  such  as  that  indicated  in  Figure  2 (b) . Two  of  these 
are  given  by  equations  (6)  and  (7) , while  the  third  is 


l 
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1 | 

J* 

*n+l  - * fa-I  ««> 

for  the  region  immediately  adjacent  to  the  mid-span,  where 
2n  is  the  mode  number.  The  negative  sign  in  equation  (24) 
is  associated  with  with  region  (n  + 1)  illustrated  in 
Figure  2(b),  while  the  positive  sign  corresponds  to  the  case 
when  the  plastic  hinge  n lies  below  the  x axis. 

Now,  substituting  equation  (24)  into  equation  (2) 
then  integrating  twice  with  respect  to  x and  satisfying  the 
boundary  conditions  M’  - Q = 0,  M = + MQ  at  x = nR+  and 
equation  (5b)  at  x = L gives 

i S- 

| nn  = L - n//2  (25) 

when  using  equations  (9)  and  (10)  for  the  plastic  hinge  n. 
Equation  (12)  remains  valid  so  that 

i 

nn  " nn-l  * 

which  when  combined  with  equation  (25)  gives 

H i = L - (1  + 1/2)  /Tt)  . 

n i i 


(26) 


(27) 


+The  + sign  is  associated  with  a plastic  hinge  n as 
illustrated  in  Figure  2 (b) , while  the  - sign  corresponds 
to  the  case  when  hinge  n lies  below  the  x axis. 
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It  is  evident  from  equations  (12) , (26)  and  (27)  that 
nn_j  = L - (j  + 1/2)/2tij»  or 

n,  = L - (n  - i + 1/2)  ^Tn  (28) 

i i 

when  n - j * i.  Thus,  equation  (14)  and  equation  (28) 
with  i = 2 predicts 

n - L/{1  + (n  - 1/2)  ✓?} 

i 

which  when  combined  with  equation  (28)  gives 

"i  * 1 + ( i - 1)  /? 

TT  1 + <n  - 1/2 ) 

where  i = 1,  2,  3,  . . . , n. 

The  response  is  now  sought  for  a fully  clamped  beam 
subjected  to  an  impulsive  velocity  field  which  has  a shape 
of  the  same  form  as  Figure  2(b)  with  equal  peak  initial 
velocities  |vQ|  at  every  hinge  1,  2,  3,  . . . , n.  In  fact, 
equations  (6)  to  (11)  remain  valid  for  antisymmetric  mode 
deformations,  except  is  now  given  by  equation  (28).  Thus, 
equations  (18)  to  (22)  remain  unchanged  and  the  maximum 
transverse  deflection  and  response  duration  are  again  given 
by  equations  (20)  and  (19),  respectively,  but  with  defined 
by  equation  (29) . 

The  foregoing  theoretical  solution  is  exact  in 


(29) 


(30) 


ll 


& 

\ 


regions  1 to  n as  discussed  earlier  for  the  symmetric  case. 

It  is  evident  for  equations  (2),  (9)  and  (24)  that  M"  has 
the  same  sign  as  $ in  the  region  nn  £ x L immediately 
adjacent  to  the  mid-span.  Thus,  the  bending  moment 
distribution  is  statically  admissible  in  this  region  because 
M'  =0  and  M = + Mq  at  x = n and  M = 0 at  x = L.  Therefore, 
the  theoretical  solution  for  an  antisymmetric  modal  response 
is  both  kinematically  and  statically  admissible  so  that  it 
is  exact  according  to  the  methods  of  plasticity  developed 
for  structures  which  undergo  infinitesimal  deflections. 
Equations  (19)  and  (20)  with  n = 1 and  rij  defined  by  equation 
(29)  agree  with  the  second  mode  solution  presented  in 
Reference  ( 1 J . 


A 


Equations  (23a)  to  (23c)  with  defined  by  equation 
(30)  again  give  the  total  angular  rotation  rates  across  the 
plastic  hinges  located  at  x = 0,  x * and  x = (i  = 2,  3, 
4,  . . . , n - 1),  respectively.  It  can  be  shown  when 
substituting  equation  (30)  into  equation  (23c)  that 

lej  = 2/3w{1  + (n  - 1/2)  /J}/L.  (31) 


(29)  , 
hinge 
(31)  which 


(26)  and 

plastic 
by  equation 
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4.  Influence  of  Finite-Displacements. 


4.1  Basic  Equations. 

It  was  shown  in  Reference  ( 3 j that 


.K 


L(p- 


pw) wdA 


i \ 

i=l  * C 


(Nw  - M) eidci 


(32) 


for  an  initially  flat  rigid  perfectly  plastic  plate  of  area 
A which  deforms  into  a number  of  rigid  regions  separated  by 
r straight  line  plastic  hinges  each  of  length  Ci*  is  the 

relative  angular  rotation  rate  across  a straight  line  hinge, 
w is  the  transverse  displacement  along  a line  hinge,  and  N 
and  M are  the  membrane  force  and  bending  moment  which  act  on 
a plane  which  is  transverse  to  the  mid-plane  of  a plate  and  passes 
through  a line  hinqe. 

Now,  it  is  straightforward  to  show  that  equation  (32) 

is  simply  an  energy  conservation  statement,  and,  in  fact,  it  has 

been  recently  rederived  by  Taya  and  Mura  ( 6 J using  a 

variational  method.  The  left  hand  side  of  equation  (32)  is 

the  external  work  rate,  while  the  term  -M8^  on  the  right 

hand  side  is  related  to  the  internal  energy  dissipation  due 

to  bending  at  the  plastic  hinge  lines,  if  w is  measured 

vertically  downwards  from  the  mid-plane  of  a plate,  then 

the  bending  moment  (M)  is  positive  at  a hinge  when  the 

material  on  the  upper  surface  of  a plate  is  stretched  and 

• • 

the  associated  value  of  0^  is  negative.  The  term  Nwd^  on 


,1 
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the  right  hand  side  of  equation  (32)  is  related  to  the 

internal  energy  dissipation  of  a plate  arising  from  the 

membrane  (in-plane)  forces  at  the  plastic  hinge  lines. 

This  can  be  demonstrated  when  using  simple  geometrical 

arguments  for  moderate  transverse  deflections  for  which  the 

in-plane  displacements  are  assumed  to  be  zero.  In  the  particular 

case  of  a rigid  region  rotating  with  an  angular  rotation  rate 
• 

6 about  a simple  support  having  w * 0,  for  example,  it 
is  evident  that  w0  is  the  axial  extension  rate  of  a hinge 
at  the  other  end. 

Thus,  the  integrand  on  the  right  hand  side  of 
equation  (32) 


Da  = (Nw  - M) 6 i (33) 

can  be  interpreted  as  the  internal  energy  dissipation 
rate  per  unit  length  of  a straight  line  hinge.  The 
explicit  form  of  the  dissipation  function  (D^)  depends 
on  the  type  of  supports  around  the  boundary  of  a plate 
and  on  the  yield  condition  for  the  material. 

If  the  maximum  normal  stress  yield  criterion 
is  selected  (see  Figure  5 of  Reference  ( 3 ] ) then  it  is 
evident  from  Figure  3 that  a membrane  force 


N 2s 

«o  T 


(34) 


I 
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is  accompanied  by  a bending  moment 


4z' 


(35) 


provided  0 £ z < H/2.  Equations  (34)  and  (35)  can  be 
combined  to  give  the  well  known  yield  condition 


M_ 


- 1 + 


(36) 


<3 


Now,  consider  a fully  clamped  span  of  length  2L 

which  remains  entirely  rigid  except  at  a central  plastic 

hinge  of  length  2c  and  at  plastic  hinges  located  at  the 

supports,  each  of  length  c.  The  axial  extension  of  the 

mid-surface  of  the  span  is  w2/L,  or  w*/2L  for  each  half  span, 

where  w is  the  transverse  deflection  at  the  central  hinge.  If  the 

axial  extension  of  each  half  of  the  span  is  divided  equally 

between  the  associated  support  and  central  hinges,  then  the 

• • 

axial  strain  rate  is  e * ww/2Lc.  Similarly,  the  curvature 
change  at  the  hinges  are  w/cL  which  gives  a curvature  rate 
k » w/cL.  Thus, 


t/ < m w/2. 


(37) 


However,  if  plane  cross-sct'wxons  remain  plane  during 

• • 4- 

deformation,  then  e • zk,  which  when  combined  with  equations 


+ * 

z is  the  location  of  the  neutral  axis  in  Figure  3.. 
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(34)  and  (37)  gives 


N/N  = w/H 
o 


which  requires  w/H  < 1 in  order  to  maintain  N/Nq  < 1. 

Consequently,  equation  (36)  is  valid  if  0 < w/H  < 1,  while 

N *=  Nq  and  M = 0 when  w/H  1.  Finally,  if  equations  (36) 

and  (38)  are  substituted  into  the  dissipation  function  (33) 

for  a beam  or  plate  with  a solid  cross-section  (N  = 4M  /H) , 

o o 

then 


K * Mq(1  + 3w2/H2)0i 


provided  w/H  < 1,  and 


D.  = 4M  W0./H 
1 o 1 


when  w/H  > 1. 


If  a beam  of  length  2L  is  subjected  to  a symmetric 
or  antisymmetric  impulsive  velocity  distribution,  then  p = 0, 
and  it  can  be  shown  for  one  half  span  of  the  beam  with  unit 
width  that  the  left  hand  side  of  equation  (32) 


fL 

- ywvi 

* n 


ywwdA  * - yLWW/3 


when  w is  defined  by  equation  (1)  and  $ is  given  in  Figures 


Equation  (39)  was  derived  for  M < 0 and  0^  > 0.  However, 
equation  (39)  if  also  obtained  for  the  other  possible  case 


when  M > 0 and  0^  < 0. 
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2(a)  and  2(b). 


4.2  Symmetric  Response  (Maximum  Normal  Stress  Yield  Criterion) 

The  total  energy  dissipation  in  one  half  of  a beam 
undergoing  a symmetric  modal  response  is 


Dt  = Dq  + Dj  + (n  - 2)Di  + Dn  , n > 2 


where  the  subscripts  refer  to  the  plastic  hinges  shown  in 
Figure  2(a)  and  is  defined  by  equation  (33).  Thus,  substituting 
equations  (17)  , (23a-d) , (33)  , (36)  and  (38)  into  equation  (42) 
gives 

°T  = 2M0W{1  + “ 1)>[1  + /2(n  - 1)  + (1  + 3/2  (n  - 1)}W2/H2]/L, 

n > 2.  (43) 

Equation  (32)  can  now  be  written  with  the  aid  of  equations 
(41)  and  (43)  in  the  non-dimensional  form 


+ 8^**  - 8,  , W,  < 1 


where 


8,  - 6 { 1 + 3/I(n  - 1)  Hi  + /2(n  - 1)  }/X 

$2  * - 6{ 1 + /2(n  - 1)}2/X 

W*  = W/H  , W*  - WH/Vq2  , t*  - tV0/H 


(47a-c) 


and  X is  defined  by  equation  (21b). 
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A theoretical  solution  of  equation  (44)  may  be 
obtained  using  the  method  of  successive  approximations. 

Thus,  a second  approximation  which  satisfies  the  initial 
conditions  W*  = 0 and  W4  = 1 at  tt  s 0 can  be  written 

W*  = t*  + &2 1*2/2  - ext*4(l  + 3B2t*/5  + 322t*2/10)/12  (48) 

provided  W*  <_  1.  The  duration  of  motion  T*  is  associated 
with  the  instant  when  W*  =*  0,  or 

1 + P2T*  - 8,T,3(1  + 382T*/4  + 3822T*2/20)/3  = 0 (49) 


and  the  corresponding  maximum  permanent  transverse  displacement 
W#f  is  given  by  equation  (48)  with  t*  = T*.  It  can  be  shown 
when  8j  = 0 that  the  theoretical  predictions  of  equations  (48) 
and  (49)  reduce  to  equations  (20)  and  (19)  , respectively, 
for  infinitesimal  displacements. 


Now,  as  remarked  previously,  the  foregoing  theoretical 

solution  is  valid  provided  < 1 in  order  to  ensure  N £ Nq 

according  to  equation  (38)  . The  case  W#  > 1 is  now  examined 

for  which  N =*  N and  M = 0. 
o 

If  equations  (40)  and  (41)  are  substituted  into 
equation  (32) , then 


£ • 

- yLWW/3  « 4M  W I ♦G./H  (50) 

° i-1  1 
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when  using  equation  (1).  It  turns  out  that  ^ > 0 when  $ = 1, 
and  0^  < 0 when  <p  = -1  (see  Figure  2 (a)),  so  that  equation  (50) 
becomes 


where 


- yLWW/3  = 4MQWeT/H 
eT  = |ej  + (n  - 2)  |0.|  + |en|. 


(51) 


(52) 


Thus,  substituting  equations  (17)  and  (23b-d)  into  equation 
(51)  gives 

W*  + ct2W#  = 0 (53) 

where 

a2  = 12(1  + /2(n  - 1)  > { 1 + 2/2  (n  - 1)}/X.  (54) 


The  initial  conditions  at  t*  = t^  associated  with  equation 

• • 

(53)  are  W*  = 1 and  W*  * where  t4i  is  given  by  equation 

(48)  with  WA  * 1 and  is  given  by  the  time  derivative  of 
equation  (48)  evaluated  at  t » t#i*  Therefore,  the  solution 
of  equation  (54)  is 


w*  “ (W#i/a)  sin(a(t*  - t#i)  } -*■  cos{a(t*  - t*^^)  } 


(55) 


which  predicts  a response  duration 


T*  * + a"‘  han”I(W*i/a) 


(56) 


i 


and  a maximum  permanent  transverse  displacement 


19 


W*f  = (1  + Wj^/a2)  1/2  (57) 

4.3  Antisymmetric  Response  (Maximum  Normal  Stress  Yield  Criterion). 

The  total  energy  dissipation  in  one  half  of  a beam 
undergoing  an  antisymmetric  modal  response  is 

dt  = D0  + Di  + <n  - DDi  r (58) 

where  the  subscripts  refer  to  the  plastic  hinges  shown  in 
Figure  2(b).  Substituting  equations  (23a),  (23b),  (30),  (31), 

(33)  , (36)  and  (38)  into  equation  (58)  gives 

°T  = MoW{1  + " 1/2)}  ^2  + /2(2n  - 1)  + {2  + 3/2  (2n  - 1)  >W2/H2J /L. 

(59) 

If  equations  (41)  and  (59)  are  substituted  into  equation  (32) 
then  equation  (44)  is  again  obtained,  except  now 

B,  = 6(1  + 3/2(n  - 1/2)  Hi  + /!( n - 1/2)  }/X  (60) 

and  B2  « - 6{ 1 + /I(n  - l/2)}2/X.  (61) 

Equation  (29)  can  be  used  to  rewrite  equations  (60)  and  (61) 
in  the  form 

B,  = 6(3  - 2nJ/(Xn*j)  (62) 

and  B2  = - 6/(Xn^*).  (63) 

It  turns  out  that  equations  (45)  and  (46)  for  the 
symmetric  case  can  be  recast  into  the  same  form  as  equations 
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(62)  and  (63)  provided  nt  is  defined  according  to  equation  (16). 
Thus,  equations  (44)  and  (47)- (4$)  are  valid  for  all  symmetric 
and  antisymmetric  modal  responses  with  Bj  and  B2  defined  by 
equations  (62)  and  (63)  and  n,  defined  by  equations  (16)  and 
(29)  for  the  symmetric  and  antisymmetric  cases,  respectively. 

It  can  be  shown  that  equation  (53)  is  also  obtained 
for  antisymmetric  modal  responses  with  ^ 1 , except  now 

a2  = 12(1  + /2(n  - 1/2)}{1  + /2(2n  - 1)}/A  (64) 

which  can  be  rewritten 

a2  = 12(2  - n )/(An  2)  (65) 

when  using  equation  (29).  However,  equation  (54)  also  takes  on 
the  form  of  equation  (65)  when  is  given  by  equation  (16) 
for  symmetric  modes.  Thus,  equations  (56)  and  (57)  are  valid 
for  all  symmetric  and  antisymmetric  modal  responses  when  a 
is  defined  by  equation  (65)  and  rij  is  given  by  equations  (16) 
and  (29)  for  symmetric  and  antisymmetric  responses,  respectively. 

4.4  Symmetric  Response  (Square  Yield  Criterion). 

It  has  been  shown  in  References  ( 2]  and  ( 3]  that 
the  theoretical  predictions  of  equation  (32)  with  a square 
yxeld  curve  relating  N and  M (i.e.,  N * Nq,  M = Mq)  provide 
a lower  bound  to  the  experimental  and  "exact"  theoretical 
maximum  permanent  transverse  displacements  of  uniformly  loaded 
fully  clamped  beams,  while  the  theoretical  predictions 
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corresponding  to  an  inscribing  square  yield  surface  (i.e., 

N = 0.618  N0,  M = 0.618  MQ)  give  an  upper  bound. 

Now,  substituting  equations  (17)  , (23a-d)  , (41)  and 
N = N0  and  M = + Mq  into  equation  (32)  gives 

W*  + <*2W*  = - a2  (1  + /2(n  - l)}/{2  + 4/I(n  - 1)}  , (66) 


where  a2  is  defined  by  equation  (54) . The  initial  conditions 
of  equation  (66)  are  W*  = 0 and  W*  = 1 at  t*  = 0.  Thus,  the 
duration  of  response  is 


1 £a“ 1 


T*  = a-1  tan-1  a-1{2  + 4/J(n  - 1)}/{1  + /2{n 


(67) 


and  the  associated  maximum  permanent  transverse  displacement  is 


W 


*f 


(1  + /2(n  - 1)}  [j 
2(1  + 2/2(n  - 1)  } j 


X{1  + 2/2  (n  - 1)  } 
3 { 1 + /2(n  - 1)  }3 


(68) 


Equations  (67)  and  (68)  with  n = 1 for  the  first  mode  respectively 
agree  with  equations  (35)  and  (34)  in  Reference  [1].  It  is 
straightforward  to  show  when  tan  aT*  = aT*  and 
X { 1 + 2/2(n  - 1 ) > { 1 + /2(n  - 1) }~3/3  <<  1 that  equations  (67) 
and  (68)  respectively  reduce  to  equations  (19)  and  (20)  for 
infinitesimal  displacements. 

4.5  Antisymmetric  Response  (Square  Yield  Criterion). 

It  can  be  shown  when  substituting  equations  (23a) , (23b) , 
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(30),  (31),  (41)  and  N = Nq,  M = + MQ  into  equation  (32) 
that 


W*  + a2W*  = 


a2 { 1 + /2(n  ~ 1/2)} 
2{ 1 + 2/I(n  - 1/2)} 


(69) 


where 

a2  = 12(1  + /2(n  - 1/2) }{1  + 2/2 (n  - 1/2) }/X  (70) 

and  w*  = 0 and  W*  = 1 when  t*  = 0.  The  duration  of  response 
is 

T = 1/cx  tan~*{2(2  - n )/a}  (71) 

* *i 

and  the  maximum  permanent  transverse  displacement  is 

w,,  = (2  - n f 1 |(l  + xn  2 (2  - n )/3}1/a  - l|/2  (72) 

I « L * 1 M J 

where  equation  (70)  can  be  rewritten 

a2  = 12(2  - n J/(\ nj)  (73) 


when  using  equation  (29) . It  turns  out  that  equations 
(67) , (68)  and  (54)  for  the  symmetric  case  can  also  be  recast 
into  the  same  form  as  equations  (71)  to  (73)  , respectively, 
provided  Hj  is  now  defined  by  equation  (16).  Thus,  equations 
(71)  and  (72)  can  be  used  for  any  symmetric  or  antisymmetric 
response  mode  provided  Hj  is  defined  by  equation  (16)  for 
symmetric  modes  (Figure  2(a))  and  by  equation  (29)  for 
antisymmetric  modes  (Figure  2(b)). 
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5.  Discussion. 

The  influence  of  material  elasticity  was  disregarded  in 
the  theoretical  solutions  in  sections  3 and  4.  This  is  a reason- 
able simplification  for  impulsive  loadings  when  the  energy  ratio 
defined  by  equation  (44)  in  Reference  [11  is  not  too  small  [see  7,8], 
Recently,  Wu  and  Witmer  [4,5]  developed  a spatial  finite-element 
and  temporal  central-difference  computer  code  JET  3C  which  can  be 
used  to  study  the  dynamic  elastic-plastic  response  of  beams.  Some 
numerical  results  for  the  first,  second  and  third  mode  cases  are 
compared  in  Figures  4 to  9 with  the  rigid-plastic  theoretical 
predictions  presented  herein  and  the  corresponding  strain-rate 
insensitive  experimental  tests  reported  in  Reference  [1]. 

It  is  evident  from  Figures  4 to  6 that  the  simplest  theoretical 
solution  which  retains  the  influence  of  finite-deflections  (equation 
(72))  gives  reasonable  agreement  with  the  numerical  elastic-plastic 
results  and  with  the  corresponding  experimental  values  except  in  the 
third  mode  case.  It  was  observed  in  References  [2]  =ind  [3]  that 
the  same  theoretical  method  which  retained  the  influence  of  geometry 
changes  gave  good  agreement  with  experimental  results  recorded  on 
uniformly  loaded  beams  and  rectangular  plates  made  from  strain-rate 
insensitive  materials.  The  numerical  values  for  the  elastic-plastic 
cases  in  Figures  4 to  6 were  estimated  from  the  deflection-time 
histories  as  indicated  in  Figures  7 and  8.  The  experimental  results 
for  the  third  mode  case  in  Figure  6 lie  below  all  the  numerical 
and  theoretical  predictions  and  is  conceivably  due  to  the  neglect 
of  rotary  inertia,  transverse  shear  effects  and  material  strain 
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hardening  which  are  likely  to  become  more  important  for  the  higher 
modal  responses.  These  effects  were  neglected  in  all  the  theoretical 
solutions  and  transverse  shear  effects  and  material  strain  hard- 
ening were  not  retained  in  the  JET  3C  calculations.  It  should  be 
noted  that  Symonds  [9]  incorporated  the  influence  of  transverse 
shear  forces  in  a yield  criterion  and  observed  that  they  were  re- 
sponsible for  a decrease  in  the  slope  at  the  mid-span  of  a beam  loaded 
dynamically.  However,  when  these  results  were  reconsidered  in 
Reference  [ LO]  , it  turned  out  that  transverse  shear  effects  caused 
the  maximum  permanent  transverse  deflections  to  increase  as  can  be 
interpreted  from  the  results  presented  in  Figure  4 of  Reference  f 10]  . 

Nevertheless,  it  is  not  known  whether  the  incorporation  of  transverse 
shear  effects  in  the  yield  criterion  would  decrease  or  increase  the 
maximum  permanent  transverse  displacements  of  the  fully  clamped  beams 
examined  in  this  article. 

The  magnitudes  of  the  non-dimensional  transverse  shear  forces 
(Q/Qq)  listed  in  Table  1 were  estimated  from  the  bending  moments 
in  the  computer  output  of  the  JET  3C  numerical  elastic-plastic 
program.  If  rotary  inertia  is  neglected,  then  Q = - dM/dx  = - 

(Mi  ~ Mi-1)  / (xi  “ xi-i) ' where  (xi  - xi_1>  is  0.25  in.  for  the 
first  and  second  mode  cases  and  0.1  in.  for  the  third  mode  case. 


+The  JET  3C  computer  program  does  have  the  capability  to  examine 
strain  hardening  effects.  However,  the  ultimate  stress  for  the 
aluminum  6061T6511  test  specimens  is  approximately  only  1.4  percent 
larger  than  the  0.1  percent  offset  yield  stress  (aQ)  fl]. 
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The  fully  plastic  transverse  shear  force  is  Q 


evident  from  Table  1 that  the  ratio  Q/Q  increases 


with  increase 


in  mode  number  despite  the  fact  the  W,/H  is  smaller  for  the  higher 


modes.  However,  in  order  to  reduce  the  amount  of  output  data  and 


lower  the  printing  costs,  the  bending  moment  distribution  across  a 


beam  was  printed  at  approximately  20  psec  intervals.  Thus,  larger 


values  of  Q/Q  could  have  occurred  at  intermediate  times!  Furthermore 


third  mode  case  in  Table  1 is  alternatively  estimated  at  x = 0.25  in 


using  the  values  of  the  bending  moments  at  x = 0.125  in.  and 


x = 0.375  in.,  then  it  turns  out  that  Q/Q  * 0.28  at  t = 20.3  psec 


which  is  still  larger  than  Q/Qq  in  the  first  and  second  mode  cases 
which  were  calculated  for  a 0.25  in.  separation  of  the  modes.  These 
observations  would  appear  to  justify  further  investigations  in  order 
to  determine  the  importance  of  transverse  shear  forces  on  plastic 
yielding  and  to  seek  the  influence  of  shear  deformations  on  the 
higher  modal  response  of  beams . 

The  axial  strains  according  to  the  JET  3C  numerical  method  are 
listed  in  Table  2 for  a few  locations,  on  the  upper  and  lower 
surfaces  of  the  beams  examined  in  Figures  7 and  8. 


+The  JET  3C  computer  program  was  used  to  repeat  the  third  mode 
case  in  Table  1 with  output  printed  at  approximately  1 psec  intervals 
for  the  first  150  pBec  of  the  response.  It  transpired  that  at 
x * 0.1  in.,  Q/Q  equalled  0.41,  0.41  and  0.40  when  t equalled  1.02 
peec,  16.25  psecand  26.4  psec,  respectively. 


The  first  and  third  mode  classical  vibration  periods  for  a 
fully  clamped  beam  with  the  same  parameters  as  the  one  considered  in 
Figure  8 are  596  usee  and  110  usee,  respectively.  These  values  are 
predicated  on  the  assumptions  that  the  response  is  entirely  elastic 
and  the  displacements  remain  infinitesimal.  It  is  evident  from 
Figure  8 that  the  period  of  the  high  frequency  vibration  of  the 
deformed  beam  is  approximately  85  usee  and  is  associated  with  a 
third  mode  vibration.  The  low  frequency  vibration  has  a period  of 
approxiamtely  595  usee  and  is  related  to  vibrations  in  the  first 
mode.  It  is  important  to  emphasize  that  the  classical  values  were 
developed  for  an  initially  straight  beam  whereas  the  numerical 
results  in  Figure  8 are  associated  with  the  vibrations  of  a per- 
manently deformed  beam. 

The  amplitude  of  the  first  modal  vibration  for  the  beam 
examined  in  Figure  8 is  sensitive  to  the  relative  magnitude  of  the 
initial  velocities  at  the  center  and  outer  hinges  (See  Figure  2(c) 
of  Reference  11]  ).  This  velocity  ratio  (VQ  at  center/VQ  at  outer 
peak)  is  1.022  for  the  beam  which  is  examined  in  Figure  8.  The 
amplitude  of  the  first  modal  vibration  is  somewhat  larger 
(w/H=0 . 07  at  the  beam  centre)  when  the  velocity  ratio  is  unity 
which  is  the  initial  velocity  distribution  required  for  a pure 
third  modal  response  according  to  a rigid-plastic  theory  based 
on  infinitesimal  displacements  (Figure  2(c)  of  [1]  ).  The  JET 
3C  computer  program  was  also  used  to  obtain  the  response  of  the  same 
beam  with  an  initial  velocity  amplitude  ratio  of  0.978  which  gave 
rise  to  the  first  modal  vibrations  with  an  amplitude  of  w/H=0.10 
at  the  beam  center. 
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The  theoretical  methods  developed  in  Reference  [1]  and  herein 
disregard  any  in-plane  or  axial  displacements  as  is  customary  for 
this  type  of  problem.  However,  it  turns  out  that  the  in-plane 
displacements  (u)  predicted  by  the  numerical  elastic-plastic  method 
are  generally  one  order  of  magnitude  smaller  than  the  associated 
lateral  or  transverse  displacements  (w)  as  shown  in  Table  3. 

The  energy  ratio  (Er>  for  a beam  is  defined  in  Reference  [1]  as 

the  ratio  between  the  initial  kinetic  energy  (K.E.)  and  the 

maximum  possible  amount  of  strain  energy  (S.E.)  which  can  be  absorbed 

by  a beam  in  a wholly  elastic  manner.  The  maximum  elastic  strain 

energy  is  estimated  for  the  present  problem  by  simply  multiplying 

the  volume  of  material  by  o e /2  where  e = o / E.  The  temporal 

1 o o o o 

variations  of  the  total  plastic  work,  total  elastic  strain  energy 
(S.E.)  and  total  kinetic  energy  (K.E.)  in  a beam  according  to  the 
JET  3C  numerical  procedure  are  shown  in  Figure  10  for  three 
different  cases.  One  suitable  energy  ratio  for  the  numerical 
results  is  given  by  the  ratio  of  the  initial  kinetic  energy  to 
the  maximum  residual  elastic  strain  energy  when  plastic  flow  ceases. 
In  this  circumstance,  the  energy  ratios  turn  out  to  be  13.3,  24.5 
and  46.6  for  the  beams  examined  in  Figures  10(a)  to  10(c),  whereas 


the  corresponding  values  of  Er  given  by  equation  (44)  in  Reference 


[1]  are  2.01,  5.24  and  5.56,  respectively.  This  confirms,  at 


least  for  these  particular  cases,  that  the  general  procedure  for 
calculating  energy  ratios  developed  in  the  Appendix  of  Reference 
(ID  and  specialized  to  the  present  problem  in  Reference  [1]  is 
conservative . 


in  twmrr  •- 
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6.  Conclusions. 

The  higher  modal  dynamic  plastic  response  of  fully  clamped 
beams  has  been  examined  using  various  rigid  perfectly  plastic 
theoretical  procedures  and  a numerical  elastic-plastic  computer 
code.  The  theoretical  predictions  of  equation  (72) , which  retains 
the  influence  of  geometry  changes,  and  the  numerical  results  agree 
reasonably  well  with  the  corresponding  experimental  values  on  strain 
rate  insensitive  beams  which  were  subjected  to  first  and  second 
modal  initial  velocity  distributions.  However,  all  the  experimental 
results  for  the  third  mode  case  lie  below  the  various  theoretical 
methods  and  numerical  results,  the  reason  for  which  is  possibly 
due  to  the  neglect  of  transverse  shear  effects. 
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+ ,o  numerical  elastic-plastic  predictions  of  JET3C 

computer  program  for  beam  described  in  the  title  of  Figure  7(a) 
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(b)  Second  mode.  Specimen  number  11  in  Table  2 of 
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(R  = 1 is  total  initial  kinetic  energy) . 


' - 


OTP 


' Tv't* 


FIGURE  4 


t ( sec  ) 
(a) 


t ( fi  sec  ) 


t ( fi  sec  ) 
(c) 

FIGURE  10 


UNCLASSIFIED 

LLUWITV  CLASSIFICATION  OF  THIS  PAOtfHhan  (Jala  Enlarad)  

Continuation  of  Box  20  on  front  side 

and  second  modal  initial  velocity  distributions.  However,  all 
the  experimental  results  for  the  third  mode  case  lie  below  the 
various  theoretical  methods  and  numerical  results,  the  reason 
for  which  is  possibly  due  to  the  neglect  of  transverse  shear  ef- 
fects . 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAOECWTlan  Data  emararf) 


